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ABSTRACT 



The gravitational lensing magnification or demagnification due to large-scale structures 



Introduction 



induces a scatter in peak magnitudes of high redshift type la supernovae (SNe la). The amplitude 
of the lensing dispersion strongly depends on that of density fluctuations characterized by the 
' cr s parameter. Therefore the value of ct§ is constrained by measuring the dispersion in the peak 

magnitudes. We examine how well SN la data will provide a constraint on the value of as using 
a likelihood analysis method. It is found that the number and quality of SN la data needed for 
placing a useful constraint on as is attainable with Next Generation Space Telescope. 
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■ Subject headings: cosmology: observations — cosmology: theory — gravitational lensing 

large-scale structure of universe — supernovae: general 

ON ' 
ON 

Qui: 

It has long been recognized that the type la supernovae (SNe la) may be a powerful tool for doing the 
cosmology for their homogeneity as well as their very high luminosity. For example, recent measurements 
of high-z SNe la have provided useful constraints on values of the cosmological parameters, present values 
of the density parameter fl m and normalized cosmological constant Q\, though they have somewhat large 
confidence interval in £l m -flA plane (Riess et al. 1998; Perlmutter et al. 1999, hereafter SCP99). There are 
now two observational projects, the Supernova Cosmology Projects^ and the High-Z Supernova Search^], 
involved in the systematic investigation of high-,2 SNe la for cosmological purposes. With their great effort, 
it seems to be promising that the quantity as well as the quality of the data are rapidly improved in near 
future. 

The SNe la are, however, not perfect standard candles but have scatter in their peak magnitudes. 
There are mainly two sources of the scatter: One is due to the intrinsic heterogeneity in SNe la which has 
been found empirically small, with a dispersion a m ~ 0.3mag in B band (Branch 1998 and references cited 
therein). Moreover it has been also pointed out that using the observed correlations between light-curve 
shape and luminosity in several different filters, the effective dispersion can be reduced to 0.12 ~ 0.17mag 
(Nugent et al. 1995; Hamuy et al. 1996; Riess, Press & Kirshner 1996). The another is the gravitational 
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lcnsing magnification (or dcmagnification) effect caused by the inhomogeneous distribution of the matter 
between SNe la and us. 

The lensing dispersion in the peak magnitudes due to large-scale structures in the cold dark matter 
(CDM) models has been investigated analytically (Frieman f997; Nakamura 1997) and numerically 
(Wambsganss et al. 1997; Wambsganss, Cen & Ostriker f998; Hamana, Martel & Futamase 1999). ft has 
been found that the dispersion depends strongly on the amplitude of fluctuations of the matter and their 
evolution, more explicitly on as, the rms fluctuation of the matter on 8ft. _1 Mpc, and on Ct m . Furthermore, 
the lensing dispersion becomes larger than O.f at redshift 0.5 if both il m and as are larger than 1. ft should 
be noted here that compact virialized objects such like individual galaxy or cluster of galaxies may also 
contribute to the lensing magnification (Holz & Wald 1998; Holz f998). However such nonlinear objects 
may cause a large magnification and the existence of lensing galaxies could be confirmed, for example, by 
deep imaging. Even if a lensing galaxy is not discovered, such an exceedingly luminous SN la should not 
be included in a normal SN la sample. Therefore, in this paper, we do not take into consideration strong 
lensing effects as a source of the scatter. 

The idea that the dispersion in peak magnitudes of the high-z SNe la may be a probe of the amplitude 
of the density fluctuations was first pointed out by Metcalf (1999). He found that the amount and quality 
of data needed for placing useful constraints on its value are attainable in a few years. This method has the 
advantage over other methods such as the two-point correlation functions of galaxies and cluster abundance 
in that the method is free from the unknown bias and the uncertain luminosity-temperature relation in 
X-ray clusters of galaxies. Metcalf (1999) parameterized the amplitude of the lensing dispersion by one 
parameter 770 which basically measures the amplitude of the appropriately projected density fluctuations 
but can not determine the values of fi m and as separately. Since as is one of most important quantities 
to study the evolution of the structures in the universe, it is worth exploring a possibility of placing a 
meaningful constraint on its value using the high-z SN la data. 

The purpose of this paper is to examine how well SN la data will place constraints on the values of 
Q m and a s in the light of the rapid increase in discovery of SNe la with redshift around or larger than 1 
in near future. For this purpose, we first re-examine the dispersion in the lensing magnifications predicted 
using, so-called, the power spectrum approach (Kaiser 1992; Nakamura 1997; Hamana et al. 1999) in §2. 
Special attention is paid to the scaling of the dispersion with as and Ct m . In §3, we investigate a possible 
constraint on cr 8 that is expected to be obtained from future SN la data using the likelihood analysis 
method, where the effect of as on the likelihood function enters through the dispersion of peak magnitudes 
due to the lensing magnifications. We also show the contour map of the likelihood function in the Vt m -as 
plane calculated using the currently available data in SCP99. Although the current data does not provide a 
useful constraint on the value of as , that will be a help to see how does two parameter degenerate in the 
plane. General discussions including the possibility of observing the SNe la at z > 1 with large telescopes 
including the Next Generation Space Telescope (NGST) are given in §4. 

2. Dispersion in gravitational lensing magnifications 

The variance of the lensing magnification, a 2 , of a point like source such like SNe la due to the 
large-scale structures can be estimated using both the Born approximation (Bernardeau, van Waerbeke 
& Mellier 1997; Schneider et al. 1998) and Limber's equation in Fourier space (Kaiser 1992; 1998). The 
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variance is related to the density power spectrum, P(k,w), by (Nakamura 1997; Hamana et al. 1999), 

Here w is the comoving radial distance, a is the scale factor defined by usual manner (e.g. Weinberg 1972) 
and normalized by its present value (i.e., ao = 1), and Jk{w) is the corresponding comoving angular 
diameter distance, defined as K^ 1 / 2 sinK~ 1 / 2 w, w, (-if) -1 / 2 smh(-K)^ 1 / 2 w for K > 0, K = 0, K < 0, 
respectively, where K is the curvature which can be expressed by K = (fl m + H\ — 1)Hq/c 2 . It should be 
noticed here that assumptions on the linear evolution and Gaussianity of the density field have not been 
used in deriving the equation ([|). We shall use the fitting formula of Peacock & Dodds (1996) (PD96 
hereafter) to describe the nonlinear evolution of density power spectra. The relationship between the 
comoving distance and the redshift z (or equivalently the scale factor a) can be derived from the Friedmann 
equation (Jain & Seljak 1997): 

w{z) = / da [n A a A + (1 - Q m - n A )a 2 + (l m a\ ~ 1/2 . (2) 

H JX/(X+z) 

In our recent paper, Hamana et al. (1999) numerically investigated the statistics of the weak 
gravitational lensing in CDM models performing the ray-tracing experiments combined with P 3 M TV-body 
simulations. We have compared the lensing dispersions obtained from the experiments with the predictions 
of the analytical approach with the PD96's fitting formula. We have found a good agreement between these 
two values within errors caused by the force resolution in P 3 M TV-body simulations. The analytical formula, 
eq. (|l]), combined with PD96's fitting formula is, therefore, a good approximation of the lensing dispersion 
for the study presented in this paper. 

We consider CDM models. The transfer function, we adopted, is given by Bardeen et al. (1986). 
Throughout this paper, we take the Hubble constant Hq = 70km/sec/Mpc which is consistent with almost 
all of the recent measurements (for a recent review, see Freedman 1999). The dispersion in peak magnitude 
of SNe la due to lensing magnification, ctgl, relates to er M by uql — 2.5 log(l + ct m ) ~ 1.0857(7^. In figure [l], 
we plot cjgl as a function of the redshift for four different sets of parameters, f2 m , Q\ and as- It is clearly 
shown in Figure 1 that ggl depends strongly on Q, rn and as but only weakly on fl\, since the effects of il\ 
on <jql enters only though the distancc-redshift relation and the growth of the power spectrum. In order 
to quantify the scaling of aQ L with as and fl m , we fit the dependence on these parameters to power laws. 
Table [l] provides such power-law fits for Einstein-de Sitter, open and J1a dominated flat cosmologies. It is 
evident from Table [l] that as has a comparable or a little weak dependence on afj L compared with that of 
fi m . Little deviation of the power of as from 2 (the relation, aQ L oc <r|, is expected for a case of the linear 
evolution of the density fluctuation spectrum) is attributed to the effect of the nonlinear evolution of the 
density fluctuations on small scales. 

3. The maximum likelihood analysis with the lensing dispersion 

Let us suppose that we observe N SNe la having a peak magnitude rrii (corrected for if-correction, 
decline rate-luminosity relation, dust extinction etc) with a magnitude error a m ,i and redshift Zj. The 
predicted magnitude-redshift relation is given by 



fK{w')f K {w - w') 
f K (w)a(w') 



dkkP(k,w') 



(1) 



m pred {z) = M + 5logV L (z,n m ,n x ), 



(3) 
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Fig. 1. — oql versus source redshift. The solid curves is for the model with (fl m , as) = (1,0,1). the 
dashed line is for (1, 0, 0.5), the long dashed line is for (0.3, 0, 1) and the long dot-dashed line is for (0.3, 0.7, 1). 
The fitting formula of Peacock & Dodds (1996) is used to describe the nonlinear evolution of the density 
power spectra. 



where M is related to the peak absolute magnitude M by 



M = M + 5 log 



and V L is the normalized luminosity distance defined by 



c/Ha 
lOpc J ' 



v L ( z ,n m ,n x ) = -±(i + z )f K (w(z)). 



(4) 



(5) 



In order to determine the parameters (in our case, fi m , Q\, M. and as), we shall maximize the Gaussian 
likelihood function defined by 



C = TT J- exp [ - 



(mi — y 



2a 



(6) 



where of = af n i + a^ L i in which we have assumed that there is no correlation between the magnitude error 
and that cased by the lensing magnification. 

We now estimate, basically following Metcalf (1999), the number of SNe la needed for a detection of the 
as value with a certain significance level. The precision with which a model parameter will be determined 
can be estimated by ensemble average of the Fisher matrix. For the case of as, that is given by 



r 2 r i / d 2 \nC \_ if [dv 2 GL ,JdasY 
-\ das 2 J 2^« i + ^ L .) 2 - 



(7) 
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Table 1 : Scaling of Oq L with erg and 0. 
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If we use the power-law fit for the lensing dispersion, i.e. <Jq L — {<t% / o^) 1 <j'q L , where quantities with the 
asterisk refer to their values in a certain model, moreover we assume that SNe la locate the same redshift and 
have the same a m , then the above equation is simplified to [0 2 8 ] _1 = N"/ag 2 aQ L /2(a 2 ri + &q L ) 2 . In Tabic 
we summarize required numbers of SNe la for 2a detection of as estimated under the above assumptions. 
The redshift is taken to be z — 0.5, 1, 1.5 and 2 with magnitude error a m = 0.15. As Table || indicates, 
it is essential to observe SNe la having a redshift larger than 1 for placing a meaningful constraint on er s . 
The above equation tells us that the number goes as a^ so it is very sensitive to this parameter. Reducing 
the magnitude error is, therefore, an another important point of this method. Similarly, N oc 0g~ 27 , thus if 
the actual value of erg is small, for example 0.5, more than ten times the number of SNe la compared with 
those in Table || is needed for the detection with the same significance level. It is, therefore, expected that 
this method will provide a strong upper limit on the value of as rather than an actual determination of the 
value. 

Table 2: Numerical values of oql and the numbers of SNe la needed for 2a detection of erg for a case of 
a s = 1. In all cases, a m is fixed to be 0.15. 
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Fig. 2. — Likelihood contours in the Vl m -as plane derived from the SN la data in Perlmutter et al. (1999). 
The contours are plotted where — 2\n£/ £ max is equal to 2.3 (solid line), 6.2 (dashed line) and 11.8 (long- 
dashed line), corresponding approximately 1, 2 and 3<r confidence contours for Gaussian likelihood function. 
The plus denotes the position of the maximum likelihood. 



Table |2j also shows that the required number is very sensitive to the cosmological model, especially to 
the density parameter, because <jgl strongly depends not only on as but also on Q m as shown in Table 
|l|. This immediately suggests that constraints obtained from SN la data will degenerate in the £l m -as 
plane. We shall investigate this point using the currently available SN la data of SCP99. We adopted 
SNe la used in 'primary fit' of SCP99 (their fit C). We also adopted the corrected peak magnitudes and 
magnitude errors summarized in Table 1 and 2 of SCP99, and thus we did not include the "stretch factor" a 
(SCP99) in the light curve-luminosity relation as a fitting parameter. The likelihood function is computed 
in four-parameter space (f2 m , SI a, M and as)- In figure ||, we plot the likelihood contours in the £l m -as 
plane, where we have not marginalized by integrating the likelihood function over other parameters 
and M) but have followed the peak, in other words, we have not used the mean but the mode. This does 
not make any significant difference as we will show below. In the lower-left region in Figure where fi m 
and as are small, no useful constraint is provided. This limitation comes from the fact that acL is smaller 
than a m for the models with a small f2 m and as- Therefore it will be the case even if we have a large, very 
high-z (z > 1) SN la sample. This limitation can be improved only by reducing a m . On the other hand, 
the upper-right region of Figure || is relatively well constrained. One may find in Figure || that the slope 
of the contour lines in the fl m -as plane are steeper than —0.5. The reason for this is that the dependence 
of @gl on VL m is stronger than that on erg as was shown in Table |l|, and the effect of f2 m on the likelihood 
functions also enters through the magnitude-redshift relation. It may be, therefore, said that SN la data 
will hardly place a lower limit on the value of as, but a future large, very high-z SN la sample can provide 
a useful upper limit on the value of as- 

One may question whether the lensing dispersion has any influence on the likelihood contours in the 
f2 m -fiA plane. In Figure 0, we plot the likelihood contours calculated with and without taking the lensing 
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Fig. 3. — Likelihood contours in the £l m -ilA plane derived from the SN la data in Perlmutter et al. (1999). 
The solid lines (dashed liens) and the plus (cross) are for the model with (without) taking the lensing 
dispersion into consideration. The contours are plotted where — 2\n£/£ max is equal to 2.3, 6.2 and 11.8. 
The cross and plus denotes the positions of the maximum likelihoods. The shaded regions are ruled out 
by other constraints: The "no big-bang" region at upper left, £l m < 0.1 is inconsistent with the amount of 
matter observed and we have simply taken Hoto > 0.6 for age. 



dispersion into account. The likelihood contours for the model without lensing dispersion are identical to 
those of SCP99 (their fit C). Figure || clearly indicates that the lensing dispersion has no significant effect 
on the constraints on n m and Q\, because the effect of the lensing dispersion on the likelihood function 
enters only through the dispersion. Therefore, the conclusion of SCP99 and also that of Riess et al. (1998) 
are not changed by the lensing dispersion due to large-scale structures. 



It may seem that the method proposed in this paper is not useful compared with, e.g., the cluster 
abundance which have provided a tight limit on the value of a combination of f2 m and a% (Eke, Cole & 
Frenk, 1996; Kitayama & Suto 1997). However one should remember that the theoretical prediction of the 
cluster abundance involves some uncertainties such like the X-ray luminosity-temperature relation and the 
bias. Our method is completely independent of the other methods in the sense that it is free from the 
relation between the distribution of dark matter and that of luminous matter, it can be a direct measure of 
us. The combined study of these methods will provide a reliable constraint in the f2 m -t78 plane. 

The most important point in using the SNe la as a probe of erg is, of course, to observe them at higher 
redshift. So far, there is no detection of SN la at z > 1. Gilliland, Nugent & Phillips (1999) detected a likely 
SN event in a revisit to Hubble Deep Field, it was associated with the galaxy at z = 1.32 (photometric), but 
no confirming spectrum of the SN was obtained. As this indicates, the main difficulty will be spectroscopy 
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of SNc. The region of spectrum that is used to do the light-curve correction redshifts to the infrared. The 
peak magnitude of a SN la at z — 1.5 is expected to be mj ~ 26 (Gilliland et al. 1999; Dahlen & Fransson 
1999). Direct spectroscopy will be very difficult for existing 8-10m telescope below 25th magnitude, but it 
will be possible with NGSlQ. A precise prediction of the number of SNe la at very high- 2 is a difficult task 
due to uncertainties in the cosmic star formation rate and the progenitor's life time. Dahlen & Fransson 
(1999) made a prediction of 75 ~ 400 SNe la per square degree down to Iab = 27 whose typical redshift 
will be z ~ 1 and have a broad redshift distribution to z ~ 2. They also predicted that 5 ~ 25 SNe la will 
be detectable per NGST field down to K' — 31.4. Therefore the number and quality of SN la data needed 
for for placing a useful constraint on as is attainable with NGST. 

We have not considered the possible evolution of SNe la properties or galactic environments which 
are of great concern for using the SNe la for cosmological purposes. If the intrinsic dispersion of the 
peak magnitude increases with redshift, the number of SNe la needed for placing a meaningful constraint 
increases rapidly. The systematic error in the peak magnitude provides incorrect constraints not only on Q m 
and Oa but also on as because these parameters are mutually related so that they have to be determined 
simultaneously. The detailed study of the possible evolutions will, of course, be a key to obtain the correct 
constraints on these parameters. The quantitative study of these issues will be done in elsewhere. 
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